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Abstract. The dynamics of wave packets in a relativistic Dirac oscillator is compared 
' to that of the Jaynes- Cummings model. The strong spin-orbit coupling of the Dirac 

C\| . oscillator produces the entanglement of the spin with the orbital motion similar to what 

is observed in the model of quantum optics. The collapses and revivals of the spin which 
" ^ ' result extend to a relativistic theory our previous findings on nonrelativistic oscillator 

fT^ , where they were known under the name of spin- orbit pendulum. There are important 

' relativistic effects (lack of periodicity, zitterbewegung, negative energy states). Many 

I of them disappear after a Foldy-Wouthuysen transformation. 
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PACS numbers: 03.65. Sq 
1. Introduction 



^ , This paper is an essay to mix together different popular models of quantum theory 



which have been developed separately and on different purposes and to present the 
Dirac oscillator as a relativistic version of the Jaynes-Cummings (JC) model with, in 
addition to the regular properties of this model, some interesting new ones related to 
the relativistic description. The Dirac oscillator, (DO) first introduced by Ito et al. 
was later shown by Cook to present unusual accidental degeneracies in its spectrum 
which were discussed from a supersymmetric viewpoint by Ui et al . It was refreshed 
later by Moshinsky and Szczepaniak and its symmetry Lie algebra was explicited by 
Quesne and Moshinsky 0. Moreno and Zentella showed also that an exact Foldy- 
Wouthuysen (FW) transformation could be performed. More recently Nogami and 
Toyama and Toyama et al have studied the behaviour of wave packets (WP) of 
the DO in the Dirac representation and in the FW representation in (1+1) dimensions. 
The aim of these authors was to study WP which could possibly be coherent. This 
reduction of the dimension was brought as an attempt to get rid of spin effects and to 
concentrate on the relativistic effects. 

Our aim is to extend the work of and to consider the full (3-1-1) dimensions 
and to show that we can there make an interesting new connection with the JC model 
as we will explain shortly. 
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The degeneracies of the eigenvalues of the DO are due to a spin orbit potential 
which is unusually large. In previous papers pj-pH^ have analyzed the time dependent 



behaviour of WP in a nonrelativistic harmonic oscillator potential with a constant spin 
orbit potential. We have shown that the behaviour of the spin shares a strong analogy 
with the observations made on the population of a two level atom in a cavity where it 
can make a two photons exchange ||12|| . 

For the relativistic DO the mechanism of collapses and revivals well known in the JC 
model is then expected to take place with some differences, [see [|l^ for a general review 



on the JC model and for a full list of references]. Since the energy spectrum is far from 
being linear the exact periodicity of the nonrelativistic oscillator is lost and the system 
evolves more on the lines of the regular JC model i. e. without exact recurrences. 
In addition to that, the spin motion should exhibit the famous zitterbewegung, this 
trembling motion should also be seen in the motion of the density of the wave as shown 
in This effect should disappear in the FW representation. 

Klein paradox is a second typical behaviour type presented by a relativistic WP 
scattered by a barrier. There it is generally said that it is a point in which Dirac equation 
cannot be interpreted as describing a single particle. This effect is more conveniently 
explained in the frame of hole theory. Studying the behaviour of a WP of the DO in 
a highly relativistic regime we have found that the WP contains a component counter- 
rotating that can indeed be interpreted with hole theory and which disappears in the 
FW representation. 

2. Summary of results on the DO 



In the following we will use the notations of fj] and of the book The time dependent 
Dirac equation is written as: 

ih = Ho \E' = c [ CK ■ (p — imujrjS) + mcj3 . (1) 
The components \&i and \E'2 of a spinor of energy E 

^2 



obey the equations 

{E^ - c^) ^1 = \c\p^ + m'u\'') - Shcumc^ - ^H^^l ■ S)] ^^ (3a) 

n 

(^2 _ ^2 ^4) ^ [^2(p2 ^ ^2^2^2) ^ 3^^^^2 ^ jm^^^ ^^^^ 

These components are thus the eigenstates of a spherical HO with a spin orbit coupling 
term respectively ±2uj/h. These large coupling strengths are responsible of the unusual 
degeneracies of the levels. The spectrum depends on a single parameter r defined as 

r = ^ . 4 
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Spectrum and degeneracies and the building of are well described in and 
For an eigenstate of energy Enij can also be labelled by n, the usual total number 
of quantas of the 3D oscillator, by the orbital and total angular momentum / and j and 
by the component m of jz. In terms of r 

Enij = mc^ VrA + 1 , (5) 
where A is defined as 

A = 2{n-j) + l if 1 = 3 -\ (6a) 

A = 2(n + j) + 3 if l = 3 + \- (Qb) 

Thus the states which obey (6a) have an infinite degeneracy. Among them those with 
n = I have the lowest value E = mc^. The states which obey [6b] have on the contrary 
a finite degeneracy. 

For the eigenvalues of eq. (3b) one should take care of the opposite sign of the spin 
orbit potential and of the different sign of the constant term. The eigenvalue is written 
in terms of n' I' and j' as 

En,i,f = mc^ VrA' + 1 , (7) 



with A' given by 



A' = 2{n-f) + 3 if I' = 3 



I 1 

2 



sa] 



A' = 2(n'+/) + 5 if r = / + |. (86) 

The need to make E and E' equal with different n and n', / and /' is fulfilled if one 
notices that \l/2 is connected to by 



I'^^) = ^ V2r(^ ■ a) |vl>i) . (9) 

E + m& 

The operator which couples \E'i and \E'2 above is expressed as the scalar product of the 
spin operator cr with a vector annihilation operator a defined by 



1 

a = ^ 



+ % 



(10) 



In order to satisfy eqs. (^^ and (|^) it is necessary that 

3=3 n' = n — 1 (11a) 

I' =1 + 1 if l = 3-\ (116) 

r =1-1 if l = 3 + \. (11c) 

These conditions imply that for the ground state of the equation satisfied by \E'i one 
has ^2 = 0. This corresponds to the fact that the DO is a supersymmetric potential 
for which and ^2 have the same spectrum apart from the absence in the spectrum 
of \I'2 of the ground state eigenvalue of 
Beside eq. (|) we have 



\^,) = I ^ V2r(o- ■ at) \^2) . (12) 

E — mc^ 
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Using = \nljm) and inserting (|^) into (|T2|) one obtains 



2^4 



[mc 



2\2 



2r \ {n'l'jm\{(T ■ a)\nljm) 



(13) 



where \n'l'jm) is a normalized harmonic oscillator state which has its quantum numbers 
defined by ( |lla| ) . From (y) and (|T3|) and using the phase conventions as in the reference 
book |T^ one obtains 



{n'l'jm\"^2) = sgn\ 



IE 



(14) 



E + mc^ 

The complex phase called sgn is defined by 

sgn= - i if /' = / + 1 = j + 1 (15a) 
sgn= +i if I' = I - I = j - I . (155) 

A normalized spinor with positive energy E = +Ep can now be expressed as 

- \nljm) 
\n'l'jm) 



2E„ 



sgn 



exp 



2En 



—iEpt 

h 



(16) 



In a similar manner a spinor with negative energy E = —Ep is written as 

- \nljm) 
\n'l'jm) 



vl/„(t) 



2Er) 



—sgn 



2 -Ed 



exp 



iEpt 



(17) 



It is interesting to note that the relative weights of the large and small components are 
formally expressed in terms of Ep exactly in the same manner as in the (1+1) model 
of [0,0. The energies are however given by eq. (|^) with conditions (6a-6b). 



3. Study of a circular wave packet. Theory 
3.1. Definition of the WP 

In the following we will study the time evolution of a circular WP of a special kind in the 
DO. For if: = we assume that the WP has a gaussian shape with an average position 
Vq and an average momentum po such that the WP moves upon a circular trajectory if 
it is left free of spin and in a nonrelativistic HO. The WP is assumed to be initially an 
eigenstate of spin with an arbitrary direction defined by two complex numbers a and f3. 
Let the normalized WP be 



^(r,0) 



(271)3/4^3/2 



exp 



2a2 



Po ■ r 





V / 



:i8) 



It is simpler to choose the axis of coordinates such that 
ro = xxo . 



(19) 
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The centroid Xq is expressed in units of the natural width of the HO and of a parameter 
called by 

xo = TiVa = ViV ^n/{muj) , (20) 
while the average momentum is taken as 

Po = ypo = yf^^/N/(T. (21) 

In (|TD|) and (^) x and y denote the appropriate unit vectors. 
The average angular momentum is then 

{L,) =xoPo = Nh. (22) 

The partial wave expansion of the WP involves only waves for which m = I and for 
which the total number of quantas of the oscillator is also /. It is given by 

/3 



l^(o)) = E^n^ = ^^^z = 

/=0 





V / 



(23) 



The weights A; are given by 



A, = (-l)'exp(-iV/2)^. (24) 



In other words (|T^) is a coherent state of the HO. We have studied in [p|-pJ]] its 
time evolution assuming that the hamiltonian is a non relativistic HO with a constant 
spin-orbit potential which depends upon a parameter called k as 

Vs.o. = K{L-a). (25) 

Associated to the DO there is a nonrelativistic hamiltonian that we will define as the 
operator in the right side of (3a) divided by 2mc^ and which thus presents the same 
degeneracies as those given by eq. (6a-b). 

3.2. Partial waves with spin up 



Let us isolate in the WP ([T8|) a partial wave / with spin up. A partial wave has two 
partners with different total angular momentum called j+ = / + | and j- = I — \ with 
respective energies and E_. is given for all the I present in (p3D by 

Ej^ = m(? in the relativistic theory (26 a) 

Ej^ = in the nonrelativistic one . (265) 

A partial wave with spin up has j = j+ and is an eigenstate of the hamiltonian with 
energy E^. It leads trivially to the phase exp{—iE^t) = exp(— ituot) in case (26a) with 
Uo = mc^/h. 
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3.3. Partial waves with spin down 



This component is coupled to the two partners and we need therefore which is given 
by 

^2 /o^^o; I i\ I 1 in the relativistic theory 

in the nonrelativistic one . 



E_ = mc^Y2r(2/ + 1) + 1 in the relativistic theory (27a) 

E_ = mc^ r{2l + 1) = huj{2l + 1) in the nonrelativistic one . (275) 

i) Let us consider first a two component spinor in the nonrelativistic theory 




V2/ + 1 

After a time t the spinor is given by 

1 / 



Vj+mj = j+ - |) + 



21 



V2/TT 



|/j_m, =j_). (28) 




^ l^^-iE+t/h _ ^-iE.t/n^ \llmi = 
(^^-iE+t/h _^ 21 e-'^-^/^) \llmi 



21 + 1 \ 

One sees that the average components of cr are given by 




At times 



t 







{ay) = 

16/ 
^ (2/ + 1) 



sin2[(2/ + l)ujt] . 



(29) 

(30 a) 
(305) 

(31) 



7r/[2cj(2/ + 1)] +mv 

the average of the spin is zero and the spin has a maximum of entanglement with the 
orbital motion. The most important component in the spinor (^9]) is for large enough 
values of / the spinor 

/ \ 



21 
21+1 



A2i+i)t ^iii^ 



(32) 



The time dependent part of the phase exp{—i2ut) in this spinor can be incorporated 
with the angular phase exp{il(j)) of the spherical harmonic. Therefore the main effect for 
this partial wave with spin down is that it rotates around Oz with angular velocity 2uj. 
Thus we have shown that a partial wave with arbitrary a and (3 contains a part which 
is at rest and a second part which rotates around Oz with angular velocity 2uj. This is 
a particular case of our previous findings |11[]: in the general case with arbitrary k the 
two waves with opposite spin move around Oz, the spin up part with angular velocity 
u — uois and the spin down one with uo + uois. In the DO we have simply uois = uj. 

ii) The corresponding relativistic spinor is written simply by adding two 
components which are zero initially. 

/ \ 

|///) 



V 








y2lTi 



/ \l3+3+ 





1)\ 



21 



J 



21 



I \IJ-J-) \ 






(33) 
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The first spinor gets the phase exp(— icjo^) while the second spinor requires an expansion 
in terms of spinors with positive and negative energy as analyzed in [0,^. The coupling 
is expressed by using eq. ([T6|) and (O) in terms of a coefficient called a; defined as 



ai 



E'i - m?c* 



2r(2/ 



\ l + 2r(2/ + l) 



'1 



where we have defined ui by 

toi = ujo ^l + 2r(2/ + l) . 
The second term of fBsI) is written at time t as 



/ \lj-J-) \ 





coscu/t — sin uit) \lj~.j-) 



It 



sgn ai sin cj/t | / — 1 j_ j_ ) / 
-^^{cosujit-i^ smuit)\lll-l) \ 



21 
21+1 



COS uJit ~ sin ujit) \lll) 



sgn ai sinuit |Z — 1 / — 1 / — 1) 




(34) 



(35) 



(36) 



(37) 



Finally reassembling the four parts of the spinor (|33|) one writes it as 
/ \ { [V2l[e-'^°^ -{cos uit-i^ sin ujit)]\l II- 1) \ 



V 



\lll) 






2/ + 1 



^iu^ot _^ 2l{cos uJit - sinujit)] \lll) 



\ 



2/(2/ + 1) sgnai sinujit\l - 1 1 - 11 - 1) 




(38) 



Comparing (^Bp to (^) we see that time evolution creates a small component in the 
relativistic spinor which is zero initially while the effect in the large components is 
through the following replacement 



-iE-t/n 



[cosuJit — I — smujit) 



-iujlt 



[1 + 



(39) 
(40) 

we obtain in the 



^1 

With similar arguments as those which lead us to eq. ([ 

relativistic case two waves which rotate in opposite sense with angular velocity 2uj. 
To get this conclusion we must linearize ui by writing 

ui = uJo[l + r{2l + 1)] = Wo + cu(2/ + 1) . (41) 

The part with positive frequency is now weighted by |(1 + ujq/uji) and the new part 
which rotates in opposite sense by ^{1 —ujq/uji). The third component of the spinor (^) 
contains also two waves with opposite sense of rotation and amplitude ai/2. Thus a 
relativistic WP with a = P = 1/V2 dissociates into three parts instead of two as in the 
nonrelativistic evolution: there is a part with spin up which does not move essentially, 
while the part with spin down is divided in two waves moving in opposite directions. 
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3.4- Spin averages 

We give below the expectation values of the operator cr for the circular WP (p^) 
assuming a = (3 = 1 / \/2 



21 + 1 
21 + 1 



1 + /(I H ^) COs(ti;/ — LjjQ)t + /(I -) CO?,{ljJi + U!o)t 

(1 + — ) sm.{ijJi - uJo)t + (1 - — ) sin{uJi + uJo)t 



2f 



(42) 
(43) 




cjf 2(2/ + l)2 
— )cosiujo-uji)t (44) 

— ) cos(co'o + ouAt 

These formulas extend to the relativistic DO those we have already discussed in 0. 
Because of the conservation of the total angular momentum there is no interference 
between the various / for the spin averages. Each partial wave depends on time because 
the energies of the spin orbit partners are different. In the nonrelativistic case the time 
factors depend only on the differences ui — ojq = {21 + l)uj. Therefore all these averages 
have period 2t^/uj. After a time called the collapse time Tc = n /{2\/2N) all the phases 
coming from all the partial waves are equally distributed and assuming high values of 
N all these averages are zero, i.e. there is a collapse of the spin! The average orbital 
angular momentum gets correspondingly an increase in order to preserve the average 
total angular momentum. This exchange was called the spin orbit pendulum since it 
occurs exactly periodically. In addition we have also shown that for a time equal to tt/uj 
the average spin is for high opposite to its initial value with the same coherence time 
around this revival. This behaviour make the spin orbit pendulum analogous to the JC 
model with the frequencies {21 + l)uj playing the role of the Rabi frequencies. 

For the DO we obtain essentially the same behaviour. However the periodicity 
is totally broken for high and even for low values of r because of the terms involving 
Uq + uji. This combination introduces in the spin motion high frequencies affecting each 
component with a different weight. The effect of this modulation is the well known 
zitterbewegung that can thus be seen in the DO on the observable of spin. Note that 
(42) and (^31) contain similar terms and similar weights, (ax) and (ay) will present then 



a similar time behaviour. The Oz component is however different since there is in (0) 
an extra term with frequency 2uJi. This term will produce two effects: an extra high 
frequency modulation and a partial revival at a time about half of the revival of the 
spin. This new revival of a purely relativistic origin concerns only the z component 
which will oscillate rapidly while the other components are zero. 
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4. Circular WP. Numerical calculations 

In order to exhibit the various effects discussed in the preceding section we have chosen 
the value iV = 20 which provides a WP well concentrated in configuration space with 
an interesting spread of its partial waves. This value was also used thoroughly in our 
previous papers on the nonrelativistic oscillator HqHTTI. For simplicity we use the units 
h = m = c = 1. Therefore our time units used in presented figures are proportional to 
uj~'^ = r~^. For nonrelativistic case the period of the motion is then T = 2n/uj = In jr. 

4-1- Spin averages in the DIRAC representation 

For the very low values of r, like r = 0.001 in Fig. |l], we are very near the nonrelativistic 
limit. One observes then a collapse of each component of the spin during an interval of 
time Tc and a revival of the spin which has lost its periodicity because of the use of the 
relativistic energies. The first revival of (ax) and (ay) occurs around time ir/uj. In the 
nonrelativistic case there was ^ a revival of (a^) with a change of sign. In the present 
case the spin rotates rapidly in the plane xOy and its maximum value is with the same 
sign as initially. The relativistic effects produce a slow decrease in the amplitude of the 
revivals. There is already a quite sensible difference in the behaviour of (a^) with time. 
This component fluctuates much more rapidly because it is richer in frequencies than 
the other two. It exhibits also a small increase at a time about half of the recurrence 
time due to these higher components. 

Each of these effects becomes more pronounced when the parameter r is given 
higher values. Figure ^ is for r = 0.025, Fig. |^ for r = 0.5. The components of the 
spin in the xOy plane oscillate much longer around each recurrence with a small period 
and the amplitudes of these recurrences decay. Again the behaviour of (cx^) is the 
most spectacular. One sees that on the average it does not get exactly to zero. The 
zitterbewegung is thus exhibited quite clearly in these time behaviour and qualitatively 
the component (cXz) differs from the other two. 

4-2. Probability densities 

We have not attempted to detect the zitterbewegung in the change of the probability 
density with time. Indeed since this effect involves high frequencies it is difficult to see 
in three dimensions. The reader is invited to read and ^ where it has been shown 
in (1+1) dimensions. The counter-rotating wave that was discussed in section p.3| is 
however easily shown for high enough values of r. In Fig. ^ the total probability density 
of the WP at the particular average radius (|19D is represented in spherical coordinates 
and for a few instants of time. What is shown was entirely explained in p.3| for each 
partial wave. A large part of the wave stays at the initial position, essentially the part 
with spin up. The wave is split into two waves which move in opposite direction and 
with the same angular velocity. They are centered around a circle with 6 = 7r/2. Two 
analyses of the WP have then been made and are represented in Fig. ^ and Fig. ^ for 
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t=10. There one finds that both of these moving parts are mainly localized in the 
second and third components of the spinor (|c2) and |c3) respectively) and that the 
counter-rotating part is almost entirely composed of negative energies. Here we are 
facing an effect totally absent from a nonrelativistic behaviour and not understood in 
a one particle theory. For lower values of r this part of the wave is hardly visible (not 
shown) . 

4.3. FOLDY-WOUTHUYSEN representation 

As derived in a FW transformation can be performed exactly on the DO. The result 
obtained is very simple and makes calculations extremely easy. The small components 
\E'2 of the Dirac representation disappears and eq. (3a) results as the only equation 
valid for The spinor has still its eigenvalues given by (|^) and the eq. (|28| ) and 
(PPP should be used for the spin averages. In other words the only relativistic effects 
are the use of these energies. The effects introduced by negative energies disappear. 
This fact was already well discussed in the (1+1) dimension model [^]. A comparison 
of the calculation in the Dirac and the FW representation enables to see exactly the 
manifestations of the zitterbewegung. Figures ^ and |^ present in FW representation 
the same cases as Figs. |l| and ^ in Dirac representation, respectively. Comparing the 
figures one see indeed that the rapid fluctuations have been washed out. In the FW 
representation the behaviour of each components of the spin is now the same. Thus it 
is the use of the relativistic energies that produce now the rapid but regular oscillations 
of the spin as well as the spread and decay of its revival. It is natural in this context to 
expect the disappearance of the component rotating in the wrong sense. The WP are 
compared in Fig. ^ at time t = 10. Only the part which rotates in the positive sense is 
left in the FW representation. 

4.4- Other spin directions 

The formulas of |3.2| and |3.3| can be combined conveniently to provide the behaviour of 
a WP pointing initially in an arbitrary direction. Such a study does not lead to a new 
dynamics. One can in this way simply put more weight on the part with spin down 
which is the most variable part. For example one can destroy almost totally the part 
not moving at the origin. Our choice of the initial direction has been made to see the 
components with spin up and down with the same magnitude. None of the other cases 
deserves a particular presentation. 

5. Summary and conclusions 

We have shown a new analogy between the relativistic DO and the JC model of quantum 
optics. The time evolution of the average spin associated to a WP in the DO is quite 
analogous to the time evolution of the occupation numbers of each of a two levels atom 
which interacts with an electromagnetic cavity. In the latter case the atom is entangled 
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with the cavity while in the former case the spin of the particle is entangled with its 
orbital angular momentum according to rules fixed by Dirac equation. In the same 
way in both models the mechanism of collapses and revivals takes place. The collapse 
of the spin is compensated by a corresponding increase of the orbital average angular 
momentum. This balance occurs periodically in the nonrelativistic case p|-p!T|. We have 
proposed the name of spin-orbit pendulum for this effect. In the relativistic case the 
periodicity is destroyed. There is then a rich behaviour of the spin components which 
are submitted to zitterbewegung. For a WP initially thrown with its velocity in the 
xOy plane the z component of the spin contains more frequencies and exhibits therefore 
most rapid oscillations. 

Related to the relativistic description we have found the presence of a counter- 
rotating wave built mainly from negative energies states. This component is particularly 
large for the geometry of the WP we have used in our paper. Such a component is well 
identified in textbooks [|1^ describing the scattering of a WP by a barrier where it 
is associated to Klein's paradox. To our knowledge it is the first time that such an 
effect has been seen in the case of WP bound in a potential. We have been able to 
observe a similar effect (however weaker) for a WP in 1+1 dimension. It is an open 
question whether this effect also exists for a WP in a Coulomb field. The smallness of 
the spin-orbit potential in this case may make observation difficult. 

Many relativistic effects are washed out in the FW representation. Our conclusions 
confirm totally those of 0,^. The coherence of the WP is also lost because of the 
nonlinear relativistic energies. Therefore coherent states of the harmonic oscillator 
generally spread. The counter-rotating wave also disappears completely with negative 
energies and the dynamics can be more interpreted with the ordinary one particle 
interpretation. This dynamics resembles then the well known dynamics of the population 
inversion of the JC model with Eberly revivals of Rabi oscillations |12]. There is an 
attempt by Toyama and Nogami to provide coherent relativistic WP of DO by using 
the inverse scattering method |16|. If these WP could be defined also for a 3+1 oscillator 
we would have then probably a relativistic spin orbit pendulum with a dynamics similar 
to the nonrelativistic one. To our knowledge these WP have not yet been constructed 
in the case where we take all dimensions into account. 
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Figure 1. Time evolution of average values of spin components for N = 20, 
r = 0.001 in Dirac representation. Note that values of (ctz) are 5 times enlarged. 
The nonrelativistic case (dashed line) is included in the upper curve for comparison 
(in this case the period is T = 27r/r = 6283.185). 



Figure 2. The same as in Fig. |i| but for r = 0.025. 



Figure 3. The same as in Fig. m but for r = 0.5. 



Figure 4. WP motion for = 20, r = 0.5 in Dirac representation. Presented is the 
total probability density j^'p on the surface of the sphere with radius r = xq. Note 
that motion of this circular WP remains close to the equator (narrow 9 range) . 



Figure 5. Contributions from positive end negative energy states for iV = 20, r = 0.5 
f = 10 in Dirac representation. 



Figure 6. Contributions from all components of the bispinor 5" (denoted as ci,C2,C3 
and C4) for TV — 20, r = 0.5 t = 10 in Dirac representation. Note different vertical 
scales and fact that the contribution from the fourth component is zero. 



Figure 7. Time evolution of average values of spin components for N — 20, r = 0.001 
in FW representation. Note that values of {(Jz) are 10 times enlarged. 



Figure 8. The same as in Fig. but for r = 0.5. 
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Figure 9. The same as in Fig. ^ but for FW representation. 



N=20 r=0.001 Dirac representation 
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N=20 r=0.025 Dirac representation 




N=20 r=0.5 Dirac representation 




N=20 r=0.5 Dirac representation 




N=20 r=0.5 t=10 Dirac representation 
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N=20 r=0.5 t=10 Dirac representation 




N=20 r=0.5 F-W representation 




